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ABSTRACT 

Context. The mechanisms giving rise to diffuse radio emission in galaxy clusters, and in particular their connection with cluster 
mergers, are still debated. 

Aims. We explore the internal dynamics of Abell 2254, which has been shown to host a very clumpy and irregular radio halo. 
Methods. Our analysis is mainly based on redshift data for 128 galaxies acquired at the Telescopio Nazionale Galileo. We combine 
galaxy velocities and positions to select 1 10 cluster galaxies and analyze its internal dynamics. We also use new (g', r' , i') photometric 
data acquired at the Isaac Newton Telescope and (V,;') photometric data available in the Subaru Archive. X-ray data from the XMM- 
Newton Science Archive are analyzed to study the hot gas component. 

Results. We estimate the cluster redshift (z) = 0.177, a high line-of-sight (LOS) velocity dispersion, crv ~ 1350 km s"' , and the 
X-ray temperature kT ~ 6.4 keV. Both our optical and X-ray analyses reveal a complex dynamical activity. The analysis of the 2D 
galaxy distribution reveals the presence of two density peaks, one at the East and the other at the West. Using the full 3D information 
we detect a high velocity (AKf Los ~ 3000 km s"'), low mass (cry ~ 200-500 km s"') group at the position of the 2D eastern 
peak. For the main system we compute a velocity dispersion cry ~ 1000-1200 km s"'. In the assumption of a bimodal system 
we estimate a mass Msys = 1.5-2.9 xlO'' /i^q Mq. The X-ray morphological analysis, based on power ratios, centroid shifts, and 
concentration parameter, confirms that Abell 2254 is a dynamically disturbed cluster. The X-ray isophotes are elongated toward the 
eastern direction, in agreement with a merger in the post core-crossing phase. A simple bimodal model finds that data are consistent 
with a bound, outgoing subcluster observed a few fractions of Gyr after the core crossing. However, both optical and X-ray analyses 
suggest that the main system is, at its time, a non relaxed structure, indicating N-S as a possible direction for a past accretion. 
Conclusions. We conclude that Abell 2254, for its mass and merging structure, fits well among typical clusters with radio halos. We 
shortly discuss as the particular irregularity of the radio halo might be linked to the complexity of the Abell 2254 structure. 

Key words. Galaxies: clusters: individual: Abell 2254 - Galaxies: clusters: general - Galaxies: kinematics and dynamics - X-rays: 
galaxies:clusters 



1. Introduction 

Merging processes constitute an essential ingredient of the evo- 
lution of galaxy clusters (see Feretti et al. 200:^ for a review). 
An interesting aspect of these phenomena is the possible con- 
nection between cluster mergers and extended, diffuse radio 
sources: halos and relics. The synchrotron radio emission of 
these sources demonstrates the existence of large-scale cluster 
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magnetic fields and of widespread relativistic particles. Cluster 
mergers have been proposed to provide the large amount of en- 
ergy necessary for electron re-acceleration to relativistic ener- 
gies and for magnetic field amplification (Tribble 1993 ; Feretti 
[19991 Feretti l20U2al Sarazin 2002 ). Radio relics ("radio gischts" 
as referred to by Kempner et al. 2004 ), which are polarized and 
elongated radio sources located in the cluster peripheral regions, 
seem to be directly associated with merger shocks (e.g., Ensslin 
et al l 1998. Roetti geret al. 1999, Ensslin & Gopal-Krishna^OOD 
Hoeft et al. l2004l l. Radio halos are unpolarized sources that per- 
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meate the cluster volume in a similar way to the X-ray emit- 
ting gas of the intracluster medium (hereafter ICM). Radio halos 
are more likely to be associated with the turbulence following a 
cluster merger, although the precise radio formation scenario re- 
mains unclear (re-acceleration vs. hadronic models e.g., Brunetti 
et al. 120091 Ensslin et al. l201 11 1. Recent semi-analytical calcula- 
tions in the framework of the turbulent re-acceleration scenario 
have allowed to derive the expectations for the statistical proper- 
ties of giant radio halos, in agreement with present observation 
that halos are found in very massive clusters (Cassano & Brunetti 
120051 Cassano et al. I2006I I. Very recently, a unified halo-relic 
model has been presented in the framework of hadronic mod- 
els where the time-dependence of the magnetic fields and of the 
cosmic ray distributions is taken into account to explain the ob- 
servational properties of both halos and (most) relics (Keshet 
1201 Oi l. In this model the ICM magnetization is triggered by a 
merger event, in part but probably not exclusively in the wake of 
merger shocks. 

Unfortunately, one has been able to study these phenom- 
ena only recently on the basis of a sufficient statistics, i.e. few 
dozen clusters hosting diffuse radio sources up to z ~ 0.5 
(e.g., Giovannini et al. 119991 see also Giovannini & Feretti 
I2OU2 I, Feretti [200 51 Ven turi et al. l2008l Bonafede et al. l2009l 
Giovannini et al. 120091) . It is expected that new radio tele- 
scopes will largely increase the statistics of diffuse sources (e.g., 
LOFAR, Cassano et al. 2010a) and allow the study of diffuse ra- 
dio emission in low X-ray luminosity clusters to discriminate 
among theories of halo formation (e.g., Cassano et al. 120051 
Ensslin et al. 1201 II) . The study of galaxy clusters with radio 
emission offers a unique tool to estimate strength and structure 
of large-scale magnetic fields and might have important cosmo- 
logical implications (see Dolag et al. l2008l and Ferrari et al. 120081 
for recent reviews). In particular, the study of clusters with radio 
halos/relics will likely contribute to quantify the effect of the 
non-thermal pressure to the estimate of mass and temperature in 
galaxy clusters (e.g., Loeb & Mao .1994; Markevitch ,2010 ) and, 
more in general, the thermal and non thermal effects of cluster 
mergers on global properties and cosmological parameters (e.g., 
Sarazin 2004). 

From the observational point of view, there is growing evi- 
dence of the connection between diffuse radio emission and clus- 
ter mergers, since up to now diffuse radio sources have been de- 
tected only in merging systems (see Cassano et al. 2010b). In 
most cases the cluster dynamical state has been derived from X- 
ray observations (Schuecker et al. 120011 Buote 120021 Cassano 
et al. i2010b) . Optical data are a powerful way to investigate the 
presence and the dynamics of cluster mergers, too (e.g., Girardi 
& Biviano'2002 ). The spatial and kinematical analysis of mem- 
ber galaxies allow us to reveal and measure the amount of sub- 
structure, and to detect and analyze possible pre-merging clumps 
or merger remnants. This optical information is really comple- 
mentary to X-ray information since galaxies and the ICM react 
on different timescales during a merger (see, e.g., the numerical 
simulations by Roettiger et al. ,1997) . In this context, we are con- 
ducting an intensive observational and data analysis program to 
study the internal dynamics of clusters with diffuse radio emis- 
sion by using member galaxies (DARC - Dynamical Analysis 
of Radio Clusters - project, see Girardi et al.i200'fl)- 

During our observational program, we have conducted an 
intensive study of the cluster Abell 2254 (hereafter A2254). 
A2254 is a very rich. X-ray luminous cluster: Abell richness 
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class = 2 (Abell et al. fT989l l: Lx(0. 1-2.4 keV)=7. 19x10'*'* h'^ 
erg s"' (Ebeling et al. 1996). Optically, the cluster is classified 
as a Rood-Sastry morphological type "B", i.e. binary (Struble & 
Rood |1987] l. The cluster redshift reported in the literature refers 
to that z = 0.178 of the brightest cluster galaxy (Crawford et 
al. fT995l l. 

Evidence for the existence of an extended, diffuse radio 
source was reported by Owen et al. (119991 ) and the presence 
of a radio halo, having a projected size of ~ 5' , is unambigu- 
ously shown by Giovannini et al. (!1999). Both the radio and the 
X-ray emission (VLA 1.4 GHz and ROSAT HRI data, respec- 
tively) show a very clumpy and irregular structure (Govoni et al. 
1200 Ibk with a radio power at 1.4 GHz of Pi.4ghz = 2.9x10^* h'^ 
W Hz"* and a large linear size of ~ 0.9 h^Q Mpc. 

We included this cluster in our DARC sample and obtained 
new spectroscopic and photometric data from the Telescopio 
Nazionale Galileo (TNG) and the Isaac Newton Telescope 
(INT), respectively. Our present analysis is based on these op- 
tical data, Subaru imaging data and XMM-Newton Science 
Archive data, too. 

This paper is organized as follows. We present optical data 
and the cluster catalog in Sect. 2. We present our results about 
the cluster structure based on optical and X-ray data in Sects. 3 
and 4, respectively. We discuss our results and present our con- 
clusions in Sect. 5. 

Unless otherwise stated, we indicate errors at the 68% con- 
fidence level (hereafter c.l.). Throughout this paper, we use 
Hq = 70 km s-i Mpc-i and hjo = //o/(70 km s"' Mpc"') in 
a flat cosmology with Q.o = 0.3 and Q.^ - 0.7. In the adopted 
cosmology, 1' corresponds to ~ 180 kpc at the cluster red- 
shift. 



2. Galaxy data and catalog 

2.1. Spectroscopic observations 

Multi-object spectroscopic observations of A2254 were carried 
out at the TNG telescope in May/October 2009 and in March 
2010. We used DOLORES/MOS with the LR-B Grism 1, yield- 
ing a dispersion of 187 A/mm. We used the 2048 x 2048 pixels 
E2V CCD, with a pixel size of 13.5 yum. In total, we observed 
4 MOS masks for a total of 142 slits. Total exposure times were 
3600 s for two masks and 5400 s for the other two masks. 

Wavelength calibration was performed using helium and 
neon-mercury lamps. Reduction of spectroscopic data was car- 
ried out using the IRAP0 package. Radial velocities were de- 
termined using the cross-correlation technique (Tonry & Davis 
1 19791 ) implemented in the RVSAO package (developed at the 
Smithsonian Astrophysical Observatory Telescope Data Center). 
Each spectrum was correlated against five templates for a variety 
of galaxy spectral types: E, SO, Sa, Sb, and Sc (Kennicutt |1992| l. 
The template producing the highest value of "R, i.e., the param- 
eter given by RVSAO and related to the signal-to-noise ratio of 
the correlation peak, was chosen. Moreover, all spectra and their 
best correlation functions were examined visually to verify the 
redshift determination. In three cases (IDs. 24, 35 and 60; see 
Table [1), we assumed the EMS AO redshift to be a reliable esti- 
mate of the redshift. Our spectroscopic catalog lists 128 galaxies 
in the field of A2254. 

- IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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The formal errors as given by the cross-correlation are 
known to be smaller than the true errors (e.g., Malumuth et al. 
LL992; Bardelli et al. [19941 ElHngson & Yee 1994; Quintana et 
al. [2000). Duplicate observations for the same galaxy allowed 
us to estimate the true intrinsic errors in data of the same quality 
taken with the same instrument (e.g. Barrena et al. 2009). Here 
we have double determinations for seven galaxies, thus we de- 
cided to apply the procedure that had already been applied in 
the above study obtaining that true errors are larger than formal 
cross-correlation errors by a factor of 2.2. For the seven galaxies 
with two redshift estimates, we used the mean of the two mea- 
surements and the corresponding errors. As for the radial veloci- 
ties estimated through EMSAO we assumed the largest between 
the formal error and 100 km s As for the whole catalog, the 
median value of the cz errors, taking into account the above cor- 
rection, is 66 km s '. 

2.2. INT observations 

Our photometric observations were carried out with the Wide 
Field Camera (WFC), mounted at the prime focus of the 2.5m 
INT telescope. We observed A2254 in g', r' and /' SDSS filters 
in photometric conditions and a seeing of ~1.6". 

The WFC consists of a four-CCD mosaic covering a 
33'x33' field of view, with only a 20% marginally vignetted 
area. We took twelve exposures of 600 s in g' filter, nine expo- 
sures of 450 s in r' filter and nine exposures of 400 s in /' filter. So 
a total of 7200 s in g' filter, 4050 s in r' and 3600 s in /' band. The 
observations were performed making a dithering pattern. This 
observing procedure allowed us to build a "supersky" frame that 
was used to correct our images for fringing patterns (Gullixson 
119921 ). In addition, the dithering helped us to clean cosmic rays 
and avoid the effects of gaps between the CCDs in the final full 
mosaic images. Another problem associated with the wide field 
frames is the distortion of the field. To match the photometry of 
several filters, a good astrometric solution is needed to take into 
account these distortions. Using the imcoords IRAF tasks and 
taking as a reference the USNO Bl.O catalog, we were able to 
find an accurate astrometric solution (rms ~0.3") across the full 
frame. Only in the North-East edge of the image, at ~17' from 
the cluster center, we could not obtain a sufficiently good astro- 
metric solution. In order to guarantee a good match between g' , 
r' and /' catalogs, this marginal region was removed from the 
master catalog. The photometric calibration was performed by 
observing the SA108 standard Landolt field (Landolt 1992 ) that 
is also calibrated in the SDSS photometric system (Smith i2002i) . 

We finally identified galaxies in our g' , r' and images 
and measured their magnitudes with the SExtractor package 
(Bertin & Arnouts 1996) and AUTOMAG procedure. Objects 
were identified by imposing that they cover a certain minimum 
area and have a number of counts above a limiting threshold tak- 
ing the sky local background as a reference. The limiting size 
and flux were 16 pixels and 1.5 standard deviation above the sky 
counts level, respectively. The selected limiting size corresponds 
to an apparent size of 1.3", which is about the minimum seeing 
size during the observations. We also performed visual inspec- 
tions of the frames in order to deal with the best combination 
of the above parameters that remove spurious objects from the 
catalogs. 

In a few cases (e.g., close couples of galaxies, galaxies very 
close to bright stars or close to defects of the CCD) the stan- 
dard SExtractor photometric procedure failed. In these cases, we 
computed magnitudes by hand. This method consisted of mea- 
suring fluxes assuming a galaxy profile of a typical elliptical 



galaxy and scaling it to the maximum observed value. The in- 
tegration of this profile provided an estimate of the magnitude. 
This method is similar to PSF photometry, but assumes an E- 
type galaxy profile, which is more appropriate in this case. 

As a final step, we estimated and corrected the Galactic ex- 
tinction Ag' ~ 0.21, A,-f ~ 0.15 and A,- ~ 0.10 using Schlegel et 
al. (1998 ) reddening and extinction maps computed from IRAS 
and COBE/DIRBE data. We estimated that our photometric sam- 
ple is complete down to g' = 23.4 (24.3), r' = 22.1 (23.3) and 
/' = 21.9 (23.3) for S /N = 5 (3) within the observed field. We 
assigned magnitudes to all galaxies of our spectroscopic catalog. 

2.3. Subaru images 

Subaru Suprime-Cam data were retrieved from the SMOKA 
archiv43- Following the same steps adopted by Nonino et al. 
( 120091 ), the images in V-band and /'-band were coadded and re- 
duced to cover a 30'x30' field of view. Seeing conditions were 
~ 0.8"- 1.0", with a total exposure time of 1200s in V-band, and 
1920s in /'-band. The astrometric solutions were obtained using 
2MASS as reference. 

Unfortunately, the observations in both filters were collected 
under non-photometric conditions, as shown from the relative 
photometric analysis. Therefore, the final coadded images were 
not calibrated using standard stars. For the use we make in this 
study, i.e. the 2D analysis of substructure, non calibrated magni- 
tudes are quite sufficient and we refer them as i'^^^ and Vobs (see 
Sect. |33|for other details). 

2.4. Galaxy catalog 

Table [T] lists the velocity catalog (see also Fig. [2]): identifica- 
tion number of each galaxy and member galaxies, ID and IDm 
(Cols. 1 and 2, respectively); right ascension and declination, a 
and 6 (J2000, Col. 3); g', r', and /' INT magnitudes (Cols. 4, 
5 and 6); heliocentric radial velocities, v = czq (Col. 7) with 
errors, Av (Col. 8). 

2.5. Individual galaxies 

No evident dominant galaxy is present in the cluster since the 
brightest galaxy in our catalog (ID. 53, r' - 16.22, hereafter 
BCG) is less than 0.5 mag brighter then the following brightest 
galaxies (IDs. 117, 102, 112 with r' = 16.64, 16.66, and 16.67 
respectively). These bright galaxies all lie in the eastern region 
of the cluster; we refer to them as BCG2, BCG3, and BCG4. 

The BCG is close to the peak of the X-ray emission and 
has a (projected) bright companion (r' - 17.01, ID. 50 here- 
after BCGc). BCGc is characterized by [Oil], H/3, and [OIII] 
emission spectral lines, strong 24-jj MlPS-Spitzer emissiorQ UV 
emissioijfl and faint radio emission (ID. 2 in Table 3 of Rizza et 
al. 2003 ). The Subaru images show that it is a spiral galaxy with 
a large bulge and an arm directed toward the BCG (see Fig. [3]). 
We find a quite large relative velocity difference with respect to 
the BCG (~ 1400 km s '). 

Figure 2 of Govoni et al. (1200 Ibl ) suggests the presence of 
a few discrete radio sources. The brightest one - with a flux of 
~ 38 mJy at 1.4 GHz - lies in the NE cluster region where the 
field is very crowded. The likely optical counterpart might be 
the cluster galaxy ID. 93 (hereafter "RG") which has an optical 

^ http://smoka.nao.ac.jp/SUPsearch 

* see at http://sha.ipac. caltech.edu/ applications/Spitzer/SHA/| 
^ http://galex.stsci.edu/GalexView/ 
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Table 1. Velocity catalog of 128 spectroscopically measured 
galaxies in the field of the cluster A2254. 



Table 1. Continued. 



ID 


IDm 




a, 6(12000) 


s' 


r' 


i' 


V 

(km 


Av 

s ) 


1 


1 


17 


17 24.62 


+ 19 41 59 6 


18.52 


17.45 


16.98 


51510 


53 


2 


2 


17 


17 25.92 


+ 1 9 39 54 7 


20.72 


19.63 


19.25 


52688 


68 


3 


- 


17 


17 25.97 


+ 19 43 95 3 


18.62 


17.27 


16.84 


72033 


64 


4 


- 


17 


1 7 96 26 


+ 19 4? 90 5 


18.90 


17.70 


17.41 


70536 


66 


5 


3 


17 


17 26.59 


+ 19 41 31 6 


18.89 


18.17 


17.96 


52741 


77 


6 


4 


17 


17 27.62 


+ 19 42 14.8 


19.31 


18.15 


17.74 


52931 


46 


7 


- 


17 


17 28.94 


+ 1 9 40 18 1 


21.59 


20.08 


19.37 


165420 


106 


8 


5 


17 


17 30.10 


+ 19 43 30 


19.88 


18.71 


18.30 


53232 


70 


9 


6 


17 


17 31.68 


+ 1 9 49 06 1 


19.28 


18.16 


17.76 


53803 


55 


10 


- 


17 


17 32.35 


+ 19 43 13 4 


20.82 


19.89 


19.46 


78076 


84 


11 


7 


17 


17 33 53 


+ 19 41 07.1 


19.68 


19.11 


18.99 


53136 


90 


12 


- 


17 


1 7 33 55 


+ 19 42 42.1 


20.26 


19.21 


18.94 


71797 


64 


13 


- 


17 


17 34.37 


+ 1 Q 49 4S 6 


17.50 


16.51 


16.12 


36622 


37 


14 


8 


17 


17 35 35 


+ 19 42 42.8 


19.70 


18.87 


18.56 


54559 


103 


15 


9 


17 


17 35.95 


+ 19 39 27.4 


20.10 


19.10 


18.81 


52006 


139 


16 


10 


17 


17 36.43 


+ 19 41 07 8 


19.71 


18.67 


18.32 


51360 


42 


17 


11 


17 


17 36.46 


+ 1 9 39 09 9 


19.32 


18.57 


18.26 


54600 


73 


18 


12 


17 


17 36.74 


+ 1 9 49 03 6 


20.80 


19.65 


19.28 


54105 


70 


19 


13 


17 


17 37.03 


+ 19 41 1 1.4 


19.44 


18.75 


18.35 


53377 


53 


20 


- 


17 


17 37.87 


+ 19 41 17.9 


21.04 


19.72 


19.36 


62831 


62 


21 


- 


17 


17 38.04 


+ 1 9 40 04 4 


20.20 


18.95 


18.72 


70432 


46 


22 


14 


17 


17 38.04 


+ 1Q 44 7 


20.37 


19.27 


18.89 


55154 


101 


23 


15 


17 


17 38.21 


+ 1 9 40 03 7 


21.76 


20.84 


19.88 


52576 


150 


24 


16 


17 


17 38.26 


+ 19 38 46 3 


21.17 


20.03 


19.47 


53964 


100 


25 


17 


17 


17 38.95 


+ 19 40 01 6 


19.37 


18.20 


17.79 


52664 


62 


26 


- 


17 


17 39.10 


+ 19 41 33 


21.20 


19.77 


19.21 


149288 


48 


27 


18 


17 


17 39.65 


+ 19 41 06 7 


20.22 


19.05 


18.79 


51500 


121 


28 


19 


17 


17 39.70, 


+ 19 38 02 8 


20.52 


19.45 


19.07 


54068 


63 


29 


20 


17 


17 39.84, 


+ 19 38 10.7 


20.67 


19.64 


19.30 


52684 


90 


30 


21 


17 


17 39.84, 


+ 19 38 10.7 


20.67 


19.64 


19.30 


54818 


90 


31 


22 


17 


17 40.03, 


+ 19 40 52.0 


19.43 


18.35 


17.98 


52821 


33 


32 


23 


17 


17 40.63, 


+ 19 37 42.6 


20.52 


19.86 


19.63 


51947 


240 


33 


24 


17 


17 40.66, 


+ 19 41 40.9 


20.33 


19.16 


18.70 


54364 


66 


34 


25 


17 


17 40.68, 


+ 19 39 52.2 


20.98 


19.84 


19.39 


54308 


103 


35 


26 


17 


17 41.02, 


+ 19 43 53.8 


20.16 


19.37 


19.06 


51559 


100 


36 


27 


17 


17 41.33, 


+ 19 39 52.2 


19.78 


19.06 


18.81 


53081 


128 


37 




17 


17 41.38, 


+ 19 41 42.0 


18.73 


18.04 


17.75 


62608 


55 


38 


28 


17 


17 41.57, 


+ 19 41 07.4 


19.01 


18.04 


17.81 


54845 


66 


39 


29 


17 


17 42.34, 


+ 19 42 39.2 


17.98 


17.18 


16.90 


54521 


77 


40 


30 


17 


17 43.01, 


+ 19 42 26.3 


19.19 


18.05 


17.58 


53870 


62 


41 


31 


17 


17 43.06, 


+ 19 41 44.2 


18.91 


17.98 


17.57 


52292 


35 


42 


32 


17 


17 43.06, 


+ 1943 51.6 


18.75 


17.69 


17.41 


53468 


41 





ID 


IDm 




Q',5(J2000) 


s' 


r' 




V 

(km 


Av 

s-') 


43 


33 


17 


17 43.20, 


+ 19 44 59.6 


19.16 


18.28 


17.99 


50532 


51 


44 


34 


17 


17 43.61, 


+ 19 39 20.5 


19.25 


18.10 


17.77 


52880 


66 


45 


35 


17 


17 43.61, 


+ 19 42 42.1 


19.58 


18.47 


18.09 


53563 


44 


46 


36 


17 


17 43.63, 


+ 19 38 37.7 


21.34 


20.46 


20.07 


51447 


99 


47 


37 


17 


17 43.70, 


+ 19 39 24.8 


20.90 


19.97 


19.65 


52533 


189 


48 


38 


17 


17 43.80, 


+ 19 41 12.1 


19.26 


18.24 


17.96 


55744 


62 


49 


39 


17 


17 43.87, 


+ 19 41 43.4 


19.42 


18.10 


17.93 


52024 


66 


50 


40 


17 


17 45.53, 


+ 19 40 40.1 


17.90 


17.01 


16.74 


54579 


53 


51 


41 


17 


17 45.62, 


+ 1941 57.1 


21.28 


20.09 


19.72 


53513 


110 


52 


42 


17 


17 45.65, 


+ 19 38 13.6 


19.31 


18.33 


17.93 


53918 


66 


53 


43 


17 


17 45.86, 


+ 19 40 48.4 


17.44 


16.22 


15.73 


52963 


41 


54 


44 


17 


17 45.86, 


+ 19 41 47.4 


19.69 


18.62 


18.23 


51347 


55 


55 


45 


17 


17 45.86, 
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spectrum typical of ellipticals. ID. 93, with an absolute magni- 
tude 21.9andaradiopowerlogPi.4GHz(WHz"') ~ 24.5, 

can be accommodated among the radio galaxies in clusters (e.g., 
Ledlow & Owen 1996; Bardelli et al. I2010I I. However, a very 
faint, extended object close to RG is revealed by the Subaru im- 
age (Fig. ID . Due to the present resolution in radio data, there is 
some doubt on the identification of RG as the (unique) respon- 
sible of the whole radio emission of the bright discrete radio 
source. Figure |4] also shows the position of a point-Uke X-ray 
source and of a MlPS-Spitzer source. 

A full discussion of individual galaxies is far from the aim 
of this study and we only give Figs. [3] and |4]in the electronic 
version of the paper 



3. Analysis of the optical data 

3.1. Member selection 

To select cluster members among the 128 galaxies with redshifts 
we perform the ID adaptive-kernel method (hereafter DEDICA, 
Pisan i [T993| and [T996l see also Fadda et al. ,1996. Girardi et al. 
1 19961 1 ■ We search for significant peaks in the velocity distribution 
at >99% C.I.. This procedure detects A2254 as a peak at z ~ 
0.177 populated by 110 galaxies considered as fiducial cluster 
members (in the range 49 237 < v < 56 757 km s"', see Fig.|5]l. 
As for the center of A2254, we adopt the position of the BCG 
[R.A. = 17'>17™45?86, Dec.=+19°40'48.4" (J2000.0)]. The 18 
non-members are 2 and 16 foreground and background galaxies, 
respectively. 

In the previous analyses of DARC clusters we generally used 
a second step in the member selection, which uses the combina- 
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04.46, 
04.61 
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tion of position and velocity information, i.e., the "shifting gap- 
per" method by Fadda et al. ( 1996 ). This procedure rejects galax- 
ies that are too far in velocity from the main body of galaxies 
within a fixed bin that shifts along the distance from the cluster 
center (adopted values are 1000 km s"' in the cluster rest-frame 
for the velocity cut and 0.6 h^^ Mpc, or large enough to include 
15 galaxies, for the spatial shifting bin). The procedure is iter- 
ated until the number of cluster members converges to a stable 
value (see Fadda et al. for other details and their Fig. 1 for an 
example). Here, the formal application of the "shifting gapper" 
procedure would reject another four galaxies just at the border 
of the distribution of galaxies in the plane of the rest-frame ve- 
locity vs. projected clustercentric distance (see the crosses in the 
Fig.|6]- top panel). These galaxies are not spatially clustered and 
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Fig. 5. Redshift galaxy distribution. The solid line histogram 
refers to the 110 galaxies assigned to A2254 according to the 
DEDICA reconstruction method. 



have the typical early-type galaxy colors. Therefore, we decide 
to not reject them from the analysis. Global cluster quantities are 
independent of this decision within the Icr error. Note that the 
"shifting gapper" procedure is completely empirical and has the 
advantage of being independent of hypotheses about the poorly 
known dynamical status of the cluster Other available proce- 
dures combining position and velocity information are generally 
based on physical assumptions about the dynamical status of the 
cluster and the consequent possibility to apply the virial mass 
and/or the knowledge of the mass profile (e.g., den Hartog & 
Katgert 1996 see Biviano et al. 2006 for a recent application). 
Therefore, the application of these procedures is not well jus- 
tified in the case of A2254 showing evidence of substructure 
(see in the following). Making an attempt in this direction, we 
use the value of the mass obtained in Sect. 15.31 Mjys, and adopt 
a one-single cluster model described by the NFW mass profile 
(Navarro et al. 1997 , see Sect. l5.3l too) to construct the upper lim- 
its to the LOS velocities at projected distance R from the cluster 
center (criterion (i) by den Hartog & Katgert |1996l l. Figure |6] 
shows as the region of allowed velocities is roughly coincident 
with that occupied by our selected member galaxies. Thus, al- 
though in the case of the substructured cluster A2254 we do not 
expect to determine the detailed cluster membership, it supports 
the consistency between our adopted galaxy membership and the 
value we obtain for the mass estimate. 



3.2. Global cluster properties 

By applying the bi weight estimator to the 110 cluster members 
(Beers et al. 1 19901 ROSTAT software), we compute a mean clus- 
ter redshift of (z) = 0.1772+ 0.0004, i.e. (v) = (53 128+125) 
km s"'. We estimate the LOS velocity dispersion, cry, by us- 
ing the biweight estimator and applying the cosmological cor- 
rection and the standard correction for velocity errors (Danese 
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Fig. 6. Top panel: rest-frame velocity vs. projected clustercentric 
distance for the 110 member galaxies of the fiducial sample 
(Fig. |5]l. Crosses indicate galaxies formally rejected as inter- 
lopers by the "shifting gapper" procedure, but here considered 
in a more conservative view. The cluster center coincides with 
the position of the BCG. Red solid (dashed) curves enclose the 
region of allowed values for member galaxies according to the 
criterion of den Hartog & Katgert (1996) and using the upper 
(lower) limit of our mass estimate computed in the following 
(see Sect. I5.3l l. Middle and bottom panels: differential (big cir- 
cles) and integral (small points) profiles of mean velocity and 
LOS velocity dispersion, respectively. For the differential pro- 
files, we plot the values for four annuli from the center of the 
cluster, each of 0.25 /i^g Mpc (large blue symbols). For the in- 
tegral profiles, the mean and dispersion at a given (projected) 
radius from the cluster-center is estimated by considering all 
galaxies within that radius - the first value computed on the five 
galaxies closest to the center. The error bands at the 68% c.l. are 
also shown. In the bottom panel, the horizontal line represents 
the X-ray temperature with the respective errors transformed in 
(Tv assuming the density-energy equipartition between ICM and 
galaxies, i.e. y8spec = 1 (see Sect.|5]l. 



et al. |1980l ). We obtain cry - 1340+^"' km s"', where errors are 
estimated through a bootstrap technique. 

To evaluate the robustness of the cry estimate, we analyze the 
velocity dispersion profile (Fig. |6ll. The integral profile shows 
an enhancement at ~ 0.4-0.5 /i^g Mpc from the cluster cen- 
ter and then flattens. This suggests that a robust value of cry is 
asymptotically reached in the external cluster regions, as found 
for most nearby clusters (e.g., Fadda et al. 119961 Girardi et al. 
I1996I I. As for the enhancement at ~ 0.4-0.5 /i^,] Mpc, it is likely 
due to the contamination of the galaxies of a secondary clump 
(see Sects. [3741 and 13. 5l i; see also, e.g., Girardi et al. d 19961 1 and 
Barrena al. (j2007) for examples in other clusters). 



Fig. 7. The 110 galaxies assigned to the cluster. Upper panel: 
Velocity distribution. The arrows indicate the velocities of the 
BCG, other brightest galaxies, and RG. Lower panel: Stripe den- 
sity plot where the arrow indicates the position of the significant 
gap- 

3.3. Velocity distribution 

We analyze the velocity distribution to search for possible devi- 
ations from Gaussianity that might provide important signatures 
of complex dynamics. For the following tests, the null hypothe- 
sis is that the velocity distribution is a single Gaussian. 

We estimate three shape estimators, i.e., the kurtosis, the 
skewness, and the scaled tail index (see, e.g.. Bird & Beers 
1993 ). There is no evidence of departures from the Gaussianity 
(see Table 2 of Bird & Beers fT993T l. 

We then investigate the presence of gaps in the velocity 
distribution. We follow the weighted gap analysis presented 
by Beers et al. ([19911 [1992', ROSTAT software). We look for 
normalized gaps larger than 2.25 since in random draws of a 
Gaussian distribution they arise at most in about 3% of the cases, 
independent of the sample size (Wainer and Schacht 1978). We 
detect one significant gap (at the 97% c.l.), which divide the 
cluster into two groups of 36 and 70 galaxies from low to high 
velocities (hereafter GVl and GV2, see Fig. [Tjl. The BCG3 is 
assigned to the GVl peak. Other bright galaxies (BCG, BCG2, 
BCG4, and BCGc) are aU assigned to the GV2 peak. 

To detect subsets in the velocity distribution we then resort 
to the Kaye's mixture model (KMM) test (Ashman et al. 1994|l. 
The KMM algorithm fits a user-specified number of Gaussian 
distributions to a dataset and assesses the improvement of that 
fit over a single Gaussian. In addition, it provides the maximum- 
likelihood estimate of the unknown n-mode Gaussians and an 
assignment of objects into groups. The KMM test is more appro- 
priate in situations where theoretical and/or empirical arguments 
indicate that a Gaussian model is reasonable. The Gaussian is 
valid in the case of cluster velocity distributions, where gravita- 
tional interactions drive the system toward a relaxed configura- 
tion with a Gaussian velocity distribution. The KMM test does 
not find a two-group partition, which provides a significantly 
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Fig. 8. Spatial distribution on the sky and relative isodensity con- 
tour map of spectroscopic cluster members, obtained with the 
2D-DEDICA method. The BCG is taken as the cluster center. 
Large square and circles indicate the location of the BCG and 
other brightest galaxies. Large triangle indicates RG. 



more accurate description of the velocity distribution than a sin- 
gle Gaussian. 

To search for a possible physical meaning of the two sub- 
clusters determined by the detected significant velocity gap, we 
also compare two by two the spatial galaxy distributions of 
GVl and GV2. We find no difference according to the the 2D 
Kolmogorov-Smirnov test (Fasano ri987l l. 

3.4. Analysis of the 2D galaxy distribution 

By applying the 2D adaptive-kernel method (2D-DEDICA) to 
the positions of A2254 galaxy members, we find one peak which 
is ~ 2' far East from the position of the BCG and the isodensity 
contours design a E-W elongated structure. (Fig. [8). 

Our spectroscopic data do not cover the entire cluster field 
and are affected by magnitude incompleteness. To overcome 
these problems, we use our photometric data samples which 
cover a larger spatial region. The INT sample has the advantage 
to have photometry available in three magnitude bands, while 
the Subaru sample allows us to extend our analysis to fainter 
galaxies. 

When more than two colors are available, it is more effective 
to select galaxies in the color-color space (Goto et al. l2002l l. In 
particular, Lu et al. ( 120091 see their Figs. 6 and 7) show as (r'- 
i') color can effectively eliminate high redshift galaxies from the 
red-sequence galaxies at z=0.2-0.3 redshift. In our INT photo- 
metric catalog we select likely members on the basis of both 
(r'-i' vs. r') and {g'-r' vs. r') color-magnitude relations (here- 
after CMRs), which indicate the early-type galaxies locus. To 
determine the CMR we apply the 2cr-clipping fitting procedure 
to the cluster members: we obtain r'-/' =0. 777-0. 021x r' and g'- 
r'=L497-0.021x r' (see Fig.|9]). Out of our photometric catalog 
we consider as likely "red" cluster members those objects lying 
within 0.1 and 0.15 mag of the CMRs {r'-i' vs. r') and (g'-r' 
vs. r'), respectively. The selected magnitude intervals are cho- 
sen with different widths to follow the different amplitudes of 
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Fig. 9. Upper panel: INT r' - i' vs. r' diagram for galaxies 
with available spectroscopy. Black circles and green crosses in- 
dicate member and no-member galaxies. For member galaxies 
large/middle/small size circles indicate late/middle/early type 
galaxies. The solid line gives the CMR determined on mem- 
ber galaxies; the dashed lines are drawn at +0.1 mag from this 
value. Lower panel: INT g' - r' vs. r' diagram for galaxies with 
available spectroscopy. The dashed lines are drawn at +0.15 mag 
from the CMR determined on member galaxies. 



the scatter in the two relations: 0.1 and 0.15 mag are round val- 
ues roughly corresponding to ~ 1.3 times the error associated 
to the fitted intercepts. Note that the difference in the scatter be- 
tween the two above CMRs is quite expected: for instance, using 
SDSS data. Goto et al. (I2002I I computed a scatter of 0.040 and 
0.081 for the two above relations, respectively. 

As a further check, we associate to each galaxy the galaxy 
type producing the highest value of H in the RVSAO proce- 
dure and divide galaxies in early-type (E, SO), middle-type (Sa, 
Sb), and late-type (Sc, Irr). The galaxies with z determined 
via EMSAO are assigned to late-type galaxies. Although these 
spectral-types do not pretend to be an alternative to a com- 
plete spectral analysis, they are useful to show how the selection 
through the CMR, in particular here the (g'-r' vs. r'), is useful 
to distinguish between early and late-type galaxies (see Fig. |9] 
- lower panel). Figure |9] shows as the selected magnitude inter- 
vals seem adequate to select early-type galaxies, thus to only use 
good tracers of the cluster substructure (e.g., Lubin et al. I2000I I 
and, above all, to avoid non-member galaxies which might bias 
our 2D analysis. 

Figure [TO] shows the contour map for the 682 likely cluster 
members having r' < 21. We find two, very significant peaks 
along the E-W direction, hereafter E- and W-peaks. The W-peak 
is at ~ r north from the BCG. A third peak lies in the WSW 
region, where our photometric data allow to study very external 
cluster regions (hereafter ext-peak). Similar results are found for 
the 780 objects with r' < 21.5. 

As for the results with r' < 21, Table |2] lists the number of 
assigned members, A^s (Col. 2); the peak position (Col. 3); the 
density (relative to the densest peak), ps (Col. 4); the value of 

for each peak, Xs (Col. 5). Ramella et al. (2007) tested the 
2D-DEDICA procedure on Monte Carlo simulations reproduc- 
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Table 2. 2D substructure from the INT and SUBARU photomet- 
ric samples. 
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Fig. 10. Spatial distribution on the sky and relative isodensity 
contour map of INT photometric cluster members with r' < 21, 
obtained with the 2D-DEDICA method. The two crosses indi- 
cate the E- and W-peaks. 



ing galaxy clusters. They show that the physical significance, i.e. 
the significance which takes into account the noise fluctuations, 
associated to the subclusters depends on the statistical signifi- 
cance of the subcluster (recovered from the value) and can 
be computed using simulations. Considering their eq. 5, the 
threshold for a sample of 682 objects is /^'fhreshoid Thus, all 
the three peaks we detect are physically significant. 

In the Subaru photometric catalog we select likely members 
on the basis of the (Vobs-'obs) 'ob.s diagram. To determine the 
CMR we apply the 2cr-clipping fitting procedure to the clus- 
ter members: we obtain yobs-'oi,s~l-910~0-01^^'obs- '^^^ 
photometric catalog, we consider as likely cluster members the 
objects lying within 0.1 mag of the CMR (see Fig.fTTTl. We stress 
that the use of the above color-magnitude relation, fully based on 
our uncalibrated magnitudes, is valid and only valid within our 
own catalog. However, we need a (rough) idea of the magnitude 
correction to understand the deepness of our Subaru samples. To 
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Fig. 11. Subaru Vohs - i'„h^ vs. /^^^ diagram for galaxies with 
available spectroscopy. Black circles and green crosses indi- 
cate member and no-member galaxies. For member galaxies 
large/middle/small size indicate late/middle/early type galaxies. 
The solid line gives the CMR determined on member galaxies; 
the dashed lines are drawn at +0.1 mag from this value. 

this aim, we compare our uncalibrated magnitudes i'^^^ and mag- 
nitudes /' of the common objects between Subaru catalog and 
the final INT catalog. The visual inspection of the i'^^^^ - 0.1 vs. 
i' plot, where 0. 1 is the Galactic extinction correction, suggests 
us that i'^^^ ~ /' + 0.6 mag. 

Figure[T2]shows the 2D-DEDICA results for the sample with 
i' < 20.5, comparable, for magnitude deepness, to the results 
from the INT sample with r' < 21. The i' < 20.5 sample con- 
firms the presence of two density peaks in the galaxy distribution 
of A2254 and the presence of the external group (841 galaxies 
in the whole field). However, while the E-peak position is close 
to the INT E-peak, the W-peak is closer to the BCG position 
than the INT W-peak and both eastern and western structures are 
someway elongated in the N-S direction. These hints for a more 
complex structure are reinforced by the analysis of the /' < 21 .5 
Subaru sample (Fig.[T2j 1257 galaxies in the whole field) where 
the W-structure is shown to be bimodal. Out of many peaks ob- 
tained through the 2D-DEDICA in the /' < 21.5 Subaru sample. 
Table |2] lists the results for the three densest peaks and the ex- 
ternal peak for comparison. Since in this case XthieshoW ~ 
we cannot formally discard the possibility of noise fluctuations, 
but Fig. 2 of Ramella et al. (i2007j) suggests a small probabil- 
ity (~ 10%) of noise contamination when the x^ is about half 
^rtireshoW show rcsults for deeper Subaru samples 

since, working on galaxies in magnitudes bins, one realizes that 
the distribution of fainter galaxies is no longer centered around 
the cluster position, a likely sign of contamination by the back- 
ground large scale structure. 

The complex structure detected by the 2D-DEDICA method 
is confirmed by the Voronoi Tessellation and Percolation (VTP) 
technique (e.g. Ramella et al. 120011 Barrena et al. 120051 1. This 
technique is non-parametric and does not smooth the data. As a 
consequence, it identifies galaxy structures irrespective of their 
shapes. For our purposes we run VTP on the same sample of 
1257 Subaru likely members with /' < 21.5. The result of the 
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Fig. 12. Spatial distribution on the sky and relative isodensity 
contour map of Subaru photometric cluster members with /' < 
20.5 and with /' < 21.5 (thin/blue and thick/red contours, re- 
spectively). Labels E and W indicate the peaks as detected by 
DEDICA in the INT sample (see Table ID. The BCG and other 
interesting galaxies are indicated as in Fig. [8] 
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Fig. 13. Galaxies belonging to structures as detected by the 
Voronoi Tessellation and Percolation technique. The algorithm 
is run on the sample of likely members with /' < 21 .5 extracted 
from the Subaru photometric catalog (see text). Open squares, 
open circles, crosses and solid circles indicate galaxies in struc- 
tures at the 95%, 98%, 99% and 99.9% c. Is., respectively. 



application of VTP is shown in Fig. [T3] VTP is run four times 
adopting four detection thresholds: galaxies identified as belong- 
ing to structures at 95%, 98%, 99% and 99.9% c.ls. are shown 
as open squares, open circles, crosses and solid circles, respec- 
tively. VTP confirms the elongation of the eastern and west- 
ern central structures in the N-S direction. In particular, the W- 
structure is clearly bimodal, and there are clues for a bimodality 
of the E-structure, too. Far from the cluster center, 17' WSW, 
VTP also confirms the existence of the ext-peak. 

Interestingly, the ext-peak can be successfully identified 
with the cluster NSC J171627-H193456 at Zphot -0.194, discov- 
ered by Gal et al. (2003, based on the galaxy catalogs from the 
digitized Second Palomar Observatory Sky Survey). However, 
the absence of spectroscopic z in the region of this system does 
not allow to infer if it is a bound companion cluster in the out- 
skirts of A2254 (projected at ~ 3 h^^ Mpc from the cluster cen- 
ter). 

3.5. 3D analysis 

3.5.1 . Significance of 3D substructure 

The existence of correlations between positions and velocities 
of cluster galaxies is a characteristic of true substructures. Here 
we use different approaches to analyze the structure of A2254 
combining position and velocity information. 

We find no evidence for a significant velocity gradient (see, 
e.g., den Hartog & Katgert 1996 and Girardi et al l 19961 for the 
details of the method). 

In order to check for the presence of substructure, we com- 
bine velocity and position information by computing the A- 
statistics devised by Dressier & Schectman (^1988 hereafter DS- 
test), which is recommended by Pinkney et al. ( 1996) as the most 
sensitive 3D test. For each galaxy, the deviation 6 is defined as 
6j = [(A^nn+l)/cry][(Vi'-V)2H-(crv,i-crv)^], where the subscript 



"1" denotes the local quantities computed over the Nan - 10 
neighbors of the galaxy. A is the sum of the 6 of the individual 
galaxies and gives the cumulative deviation of the local kine- 
matical parameters (mean velocity and velocity dispersion) from 
the global cluster parameters. The significance of A, i.e. of sub- 
structure, is checked by running 1000 Monte Carlo simulations, 
randomly shuffling the galaxy velocities. We find no significant 
presence of substructure. 

Following Pinkney et al. ( 119961 see also Ferrari et al. l2003l l. 
we apply two more classical 3D tests: the e-test (Bird, I1993I I 
based on the projected mass estimator and the centroid shift or 
ff-test (West & Bothun il990l ). The details of these tests can be 
found in the above papers. We only point out that we consider 
ten as the number of the nearest neighbors for each galaxy and 
we use the above Monte Carlo simulations to compute the sub- 
structure significance. In both cases we do not find evidence for 
significant substructure. 

We also consider two kinematical estimators alternative to 
the S parameter of the DS-test, i.e. we consider separately the 
contributes of the local mean - [(Mm + 1)/o"y](Vi- V)^], and 
dispersion 6^ - [(Nnn + l)/o"v ](o"v,i - ctv)^] (see, e.g. Girardi 
et al. 119971 Ferrari et al. 120031 ). When considering the 6s esti- 
mator, we find evidence for peculiar local velocity dispersion at 
the 97% C.I.. Figure [14]- lower panel - shows the distribution 
on the sky of all galaxies, each marked by a circle: the larger the 
circle, the larger the deviation 6sj of the local velocity dispersion 
from the global cluster value. The two largest deviations are ob- 
tained for the subgroups associated to two galaxies just south of 
the BCG for which crv.i ~ 400 km s Other deviations interest 
the eastern region characterized by large values of the velocity 
dispersion. To better investigate we repeat the DS-test increasing 
the number of neighbors. As for the 6^ indicator, we obtain large 
c.l. values for peculiarity up to a 99.99% c.l. in the A^nn =25- 
50 range (see Fig. [15]- upper panel). This result suggests that 
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Fig. 14. Spatial distribution of the 110 cluster members, each 
marked by a circle. The plot is centered on the cluster center 
Upper panel: the cluster velocity field: the larger the circle, the 
larger is the galaxy velocity. Lower panel: the result of the (mod- 
ified) DS-test: the larger the circle, the larger is the deviation <5s,, 
of the local velocity dispersion from the global velocity disper- 
sion. Thin/blue and thick/red circles show where the local ve- 
locity dispersion is smaller or larger than the global value. The 
cross indicates the position of the E-peak as detected in the 2D 
analysis (INT sample). 



the peculiarity interests a large part of the cluster To individuate 
how many galaxies are involved in the substructure, we resort to 
the technique developed by Biviano et al. (120021) . who used the 
individual 5,-values of the DS method. The comparison of the 
(5s,;-values of all 1000 Monte Carlo simulations and the observed 
values shows how many galaxies are involved in the structures 
with peculiar kinematics (Fig. [15]- lower panel). 

3.5.2. Identifying subclusters in line eastern region 

Comparing the two subsamples obtained by a rough division of 
A2254 (e.g., the eastern and the western subsamples with re- 
spect to the position of the BCG), we find that the eastern sam- 
ple is characterized by a larger value of cry (see Table [3] and 
Fig. [16]- upper panel). The F-test confirms that the velocity dis- 
persions of the two samples differ at the 99.4% c.l.; no difference 
is shown between the mean velocities according to the means- 
test (see,e.g.. Press et al. 1992 for these tests). 

Figure [15]- lower panel shows as the region with high o"v,i 
values is centered roughly around the E-peak as determined in 



Fig. 15. Upper panel: the same of Fig. [14]- lower panel but for 
the DS-test with a larger number of neighbors (A^nn = 40). The 
large, faint circle encloses the region within 1.5' from the E- 
peak: galaxies inside defines the E1.5-sample. Lower panel: the 
distribution of dsj deviations of the above DS-test. The solid line 
represents the observations, the dashed line the distribution for 
the galaxies of the 1000 simulated clusters, normalized to the 
observed number. 



the 2D analysis. We select the galaxies within a circular region 
centered in the position of the E-peak. Using a radius of 1.5', 
comparable to the distance between E- and W-peaks, we obtain 
a subsample of 27 galaxies, the E1.5-sample, characterized by 
a large value of cry ~ 1700 km s"'. When applying the ID- 
DEDICA method to this subsample we detect two significant 
peaks with 21 and 6 galaxies at ~ 52 500 and 56 300 km s"', 
respectively (see E1.5LV and E1.5HV in Table[3]l. Figure [T6l- 
lower panel - shows the results of the ID-DEDICA method for 
the E1.5-sample, as well as for other two samples defined by a 
larger and a smaller radius (2' and 1'). The choice of a large 
radius gives one-(asymmetric) peak in the galaxy density distri- 
bution in the velocity space, while a too small radius gives a list 
of non-significant peaks due to the poor statistics involved. 

3.5.3. Identifying subclusters through the Htree method 

We resort to the method devised by Serna & Gerbal ( 119961 here- 
after Htree-method; see, e.g., Durret et al. 2010 for a recent ap- 
plication). This method uses a hierarchical clustering analysis to 
determine the relationship between galaxies according to their 
relative binding energies. The method assumes a constant value 
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Table 3. Kinematical properties of the whole system and galaxy subsystems. 
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Fig. 16. Upper panel: velocity distributions of the galaxies in the 
western and eastern samples (thick/red and thin/black lines, re- 
spectively). The velocity distribution of a subsample of the east- 
ern sample, i.e. the 27 galaxies within a radius of 1.5' from the 
E-peak (E1.5-sample), is also shown (black thick line). Lower 
panel: the galaxy density distribution in the velocity space, as 
provided by the DEDICA adaptive-kernel reconstruction method 
for the E1.5-sample (solid line) and for two samples contained 
within a larger and a smaller radii, E2.0- and El.O-samples 
(dashed and dotted lines, respectively). Unit on the y-axis is nor- 
malized to the density of the highest peak. 



for the mass-to-light ratio of galaxies and Serna & Gerbal sug- 
gest a value comparable to that of clusters. Here we take a value 
of M/Lr=150 /i7o Mq/Lq as suggested by large statistical stud- 



ies (e.g., Girardi et al. 12000 ', Popesso et al. 120051) . The (gross) 
results are quite robust against the choice of the value of M/Lr- 

Figure[T7]shows the resulting dendogram, where the total en- 
ergy appears horizontally. Galaxy pairs and subgroups of galax- 
ies then appear with a lower total energy. Starting from the less 
deep part of the energy levels (right part of the dendogram) we 
find HTl, which we interpret as the main system, and the group 
HT2. HT2 is a high-velocity eastern group and contains all mem- 
bers of E1.5HV. HTl is a "substructured" main system and is 
formed by HTll, which has the BCG in its potential well, and 
HT12, i.e. a few low velocity galaxies located in the western re- 
gion (in particular, the BCG3). Note, however, that at the bottom 
of the binding energy of HTl 1 we find HTl 111, which is not a 
core in the classical sense. In fact, HTllll contains HTl 11 11 
and HTl 1 1 12. HTl 1 1 12 is formed by three galaxies close to the 
cluster center; the BCG and BCGc pair, and a small compact 
galaxy located between them. HTlllll is close to the E-peak 
and is formed by four galaxies, among which the BCG2 and 
BCG4. 



3.5.4. Identifying subclusters tlirougli line 3D-DEDICA 
metliod 



Finally, we apply the 3D-DED1CA method (Bardelli et al. [T998] l. 
Table [3] lists the five most significant (^ 94% c.l.) groups 
(DED1-DED5) composed by at least three galaxies. The picture 
revealed by the 3D-DED1CA test is that of a very complex struc- 
ture, with three and two subclusters in the eastern and western 
regions, respectively. 



Figure [TSl- left col. - shows the spatial distribution of the 
groups detected by above methods and Table[3]lists the kinemati- 
cal properties: the number of assigned members, A^g (Col. 2); the 
mean velocity, (v> (Col. 3); the velocity dispersion with boot- 
strap eiTors, cry (Col. 4); a short comment summarizing our 
interpretation (Col. 5). Following properties refer to the corre- 
sponding KMM subclusters, see below (Cols. 6, 7, and 8). 
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V (10^ km s-') <52 



Fig. 17. Dendogram obtained through the Serna & Gerbal d 19961 1 
algorithm. The abscissa is the binding energy (here in arbitrary 
unit with the deepest negative energy levels on the left) while 
the catalog numbers of the various member galaxies are shown 
along the ordinate (IDm in Table [T]i. 



3.5.5. 3D-KMM method 

We use the above results to apply the full 3D-KMM method, 
which, starting from the above seeds, allow us to assign all 
galaxies to subclusters. 

Using E1.5LV and E1.5HV as seeds for a initial two-group 
partition, we find that a partition of 103 and 7 galaxies is a signif- 
icantly more accurate description of the 3D galaxy distribution 
than a single 3D Gaussian (at the 99.7% c.l.). We note that the 
KMM-E1.5LV subcluster i.e. the main system, has a significant 
weighted gap in its velocity distribution (as in the whole sample, 
see Sect. 13.31 ) but does not show any trace of peculiar kinemat- 
ics according to the DS-test. Using HTl, HT2, and HTllll as 
initial seeds, the 3D-KMM converges to a three-group partition 
significant at the 97.5% c.l.. Finally, using DED1-DED5 groups 
as initial seeds we find the 3D-KMM test converges to a five- 
group partition significant at the 99.9% c.l.. 

Figure [18]- right col. - shows the spatial distribution of the 
obtained partitions. 



4. X-ray morphological and spectral analysis 

4.1. Observation and data reduction 

A2254 was observed by XMM-Newton on August 26, 2009. The 
data were reduced with SAS vlO.O using the tasks emchain and 
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Fig. 18. Spatial distribution on the sky of the galaxies of the sub- 
clusters detected through our 3D analysis and the cluster par- 
titions obtained through the 3D-KMM method using them as 
seeds {left and right panels, respectively). In the left panels sym- 
bols refer to different groups. The different colors refer to differ- 
ent mean velocities according to the legend at the top of the fig- 
ure. In the right panels symbols indicate subclusters found using 
as seeds the respective groups on the left and colors refer again 
to the velocity according to the above legend. 



epchain. We only considered event patterns 0-12 for MOS and 
for pn, and the data were cleaned using the standard proce- 
dures for bright pixels and hot columns removal (by applying 
the expression FLAG -- 0) and pn out-of-time correction. The 
energy scale of the pn over the whole spectral bandpass was 
further improved with the task epreject. Periods of high back- 
ground event counts due to soft protons were filtered using a 
threshold cut method: we performed a Gaussian fit to the peak 
of the histogram of the light curve and excluded the time bins 
where the count rate lies more than 3cr away from the mean. 
The light curves were extracted from regions of least source 
contamination (excising the bright object core in the central 8' 
and the point source list from the SOC pipeline, after visual in- 
spection) in two different energy bands: a hard band, 10-12 keV 
for MOS and 10-13 keV for pn (using 100s bins), and a wider 
band, 0.5-10 keV (using 10s bins), as a safety check for possi- 
ble flares with soft spectra. The flaring periods thus determined 
were further checked by visual inspection of the light curves, 
resulting in net exposures of 57.0 ks, 56.8 ks, and 45.1 ks for 
the MOSl, MOS2, and pn detectors, respectively. Point sources 
were detected using the task ewavelet in the energy band 0.5-10 
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keV and checked by eye on images generated for each detector. 
Detected point sources from all detectors were merged, and the 
events in the corresponding regions removed from the event list 
using circular regions of 25" radius centered at the source po- 
sition. The area lost by point source exclusion, CCD gaps, and 
bad pixels was calculated using a mask image. Redistribution 
matrix files (RMFs) and ancillary response files (ARFs) were 
generated with the SAS tasks nnfgen and arfgen, the latter in ex- 
tended source mode. Appropriate flux-weighting was performed 
for RMFs, and for ARFs using exposure-corrected images of the 
source as detector maps (with pixel size of 1 ', the minimum scale 
modeled by arfgen) to sample the variation in emission, follow- 
ing the prescription of Saxton & Siddiqui (2002). For the back- 
ground subtraction, we adopted a complete modeling of the vari- 
ous background components as described in detail in Gastaldello 
et al. ( 120071 ). We accurately modeled the actual sky background 
using a region free of cluster emission at r > 9', where the 
surface brightness profile in the 0.7-1.2 keV reaches the back- 
ground level. The spectra of the out-of-field-of-view events of 
CCD 4 and 5 of MOS 1 showed an anomalously high flux in the 
soft band (see Kuntz & Snowden 2008 ) so they were excluded 
from our surface brightness and spectral analysis. 



4.2. X-ray image and surface brightness analysis 

For each MOS detector, we create images in the 0.5-2 keV band 
and exposure maps at the monochromatic energy of 1 keV and 
then combine the images into a single exposure-corrected im- 
age, smoothed on a scale of 12", as shown in Fig[T9] We process 
the image to remove the point sources using the CIAO tool dm- 
filth, which replaces photons in the vicinity of each point source 
with a locally estimated background. We consider only the re- 
gions exposed by more than 65% of the maximum exposure (see 
Fig. |20] | in order to perform the X-ray morphology analysis. We 
adopt quantitative measures by applying the power ratios tech- 
nique (Buote & Tsai l 19961 ), the centroid shifts (Mohi- et al. |1993l l 
and the concentration parameter (Santos et al. l2008i) . 

Power ratios derive from multipole expansions of the two- 
dimensional projected gravitational potential and they can be ap- 
plied directly on X-ray images since under normal circumstances 
the X-ray emissivity can be used as a proxy for the gravitational 
potential. Power ratios are computed within a fixed circular aper- 
ture and they can be computed on the centroid of the cluster or 
on the X-ray peak. In the former case Pi/Po is null, P2/P0 gives 
information on the ellipticity of the cluster, P3 /Pq is related to a 
bimodal structure and P4/P0 is sensitive on small scale substruc- 
ture. In the latter case Pj'^/Pjj'' is non null if the cluster does not 
exhibit reflection symmetry about two orthogonal axes originat- 
ing on the peak of the emission and it is essentially a circularly 
average centroid shift. We calculate the power ratios within an 
aperture of 500 kpc, given the high signal-to-noise within this 
radius and for comparison with previous work. 

Centroid shifts indicate that the center of mass of the X-ray 
emitting gas varies with radius. Centroid shifts and power ra- 
tios are both capable of identifying highly disturbed systems or 
systems with significant, well defined substructures (Poole et al. 
120061 ). However the centroid shifts parameter, w, is more sensi- 
tive to subtler disturbances. Following the method of O'Hara et 
al. (12006 , see also Poole et al. 2006), the centroid shift is com- 
puted in a series of circular apertures centered on the cluster X- 
ray peak. The radius of the apertures is decreased in steps of 5% 
from 500 kpc to 25 kpc. 




Fig. 19. Mosaic of the MOSl and MOS2 images in flie 0.5-2 
keV energy band smoothed on a 12"scale. The image was di- 
vided by the summed exposure maps to correct for exposure 
variations. The detected point sources within the central regions 
of the cluster discussed in the text are highlighted by black cir- 
cles. 




Fig. 20. X-ray image of the diffuse emission of the A2254 clus- 
ter where point sources were removed and replaced by a local 
estimate of the surrounding emission and only regions exposed 
by more than 65% of the maximum exposure is considered. The 
contour levels of the radio VLA image (Govoni et al. l2001b] l are 
superimposed, too 



The concentration parameter, c, is defined as the ratio of 
the peak (calculated within 100 kpc) over the ambient surface 
brightness (calculated within 500 kpc). We apply these tech- 
niques to the image reported in Fig. |20] and the results are re- 
ported in Table |4] 

Following Ettori et al. (2010) we extract the surface bright- 
ness profile in the energy band 0.7-1.2 keV, in order to keep the 
background as low as possible with respect to the source. For 
this reason, we avoid the intense fluorescent instrumental lines 
of Al (~ 1.5 keV) and Si (~ 1.8 keV). We account for the X- 
ray background by including a constant background component. 
The data are grouped to have at least 20 counts per bin in order 
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Table 4. Results of the morphological analysis of the X-ray im- 
age of A2254 using power ratios, centroid shifts and concentra- 
tion. 



Pf/Po P2/P0 P3/P0 P4/P0 w c 

(xlQ-^) (xlQ-^) (xlQ-^) (xlQ-'*) (500h7„'kpc) 

1.84 + 0.47 4.25 + 1.26 2.51 +0.82 1.29+ 1.12 0.037 + 0.001 0.08 




0.1 1 

Radius (arcmin) 



Fig. 21. Surface brightness profile of the X-ray emission of 
A2254. Data from MOSl and MOS2 are plotted in black and 
red, respectively. The best fit yS-model and ratio of data over the 
model are also shown. The background level is shown as the 
black dashed line in the lower panel. 

to apply the;t^^ statistics. The fitted model is convolved with the 
XMM-Newton PSF. The joint best-fit yS-model has a core radius 
of rc = (507 + 36) h'^ kpc (169"+12") and yS = 0.95 + 0.06 for 
a;(f2/d.o.f. = 215/179 (see FigEB- 

4.3. Spectral analysis 

We extract spectra for each detector in 5 concentric circular an- 
nuli located at the the peak of the X-ray emission with bounding 
radii 0'-0.5', 0.5'-1.5', 1.5'-2.5', 2.5'-3.5', 3.5'-4.5'. The spec- 
tra from the three detectors are re-binned to ensure a signal-to- 
noise ratio of at least 3 and a minimum 20 counts per bin (nec- 
essary for the validity of the minimization method). We fit 
the spectra with an APEC thermal plasma modified by Galactic 
absorption fixed at 4.7 x 10^" cm"^ (Kalberla et al. l2005l l to each 
annulus. A correction file derived from the background mod- 
elhng is used in XSPEC as in Gastaldello et al. (.2007 J . The 
free parameters are temperature, normalization (proportional to 
emission measure), and the metallicity, with the solar abun- 
dance units of Grevesse & Sauval (1998). The spectra are ad- 
equately fitted by the model, with reduced ranging from 1.06 
in the inner annulus to 1.5 in the outermost annulus. The az- 
imuthally averaged, projected temperature and abundance pro- 
file thus obtained are shown in Fig. |22] The temperature pro- 
file is rather flat and in particular there is no evidence of a de- 
crease in the core, as found in non-cool core clusters. We cal- 
culate the mean temperature of the cluster as the temperature 
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Fig. 22. Profiles of the azimuthally averaged, projected temper- 
ature and metal abundances of A2254. The error bars are at the 
68% C.I.. The horizontal line shows the value of kTouT, i e. the 
value of kT estimated within the 0.05/? 1 go < r < 0.2/? 1 go region, 
the two radii being indicated by the two vertical lines. 

obtained with a spectral fit in the region 0.05/?igo < r < 0.2/?igo, 
where /?i8o = 1780(fer/5keV)'/^/i(z)"' h'^ kpc (Arnaud et al. 
2005 . with h{z) - (Qm(1 + zf + ^a)'/^) using an it erative pro- 
cedure to calculate kT and R^q (Rossetti & Molendi llOTOb . We 
find kT = 6.38 + 0.25 keV and /?igo = 2.08 h'^^ Mpc. The value 
of kT thus estimated within the 0.05/?i8o < r < 0.2/?i8o region, 
hereafter kTouT, is a good proxy for the global temperature (see 
Fig. 4 of Leccardi & Molendi 2008 where kTy[ is the global tem- 
perature). In particular as for A2254, Figure |22] shows as the 
measure of kT has already large uncertainties around 0.4/?i8o- 

To characterize temperature and emission measure gradients, 
we compare the central to a global value of the pseudo-entropy 
measure. Following Leccardi et al. (1201 Oi l, we compute the 
pseudo-entropy ratio £ - (kTn^ / kTouj) x (£Min/-£Mout) 
where kT and EM are the X-ray temperature and emission 
measure, and IN and OUT define the region within the radius 
~0.05 Rm) and the region within the the annulus with bounding 
radii 0.05-0.20 /?i8o, respectively. Note that Rossetti & Molendi 
(2010) proved that the pseudo-entropy ratio used in their study 
is a good proxy of the deprojected quantity and Rossetti et al. 
( 201 1. their Fig. 3) found a good correlation between the pseudo- 
entropy ratio and, for example, the central entropy defined in 
Cavagnolo et al. (2009). We measure £ = 0.833 ± 0.097 which 
places A2254 among clusters with high entropy core (see also 
Sect. 1511). 

We also prepare a two dimensional map of X-ray temper- 
ature, starting from the EPIC images. We use a modified ver- 
sion of the adaptive binning and broad band fitting technique 
described in Rossetti et al. ( 120071 ). where we substitute the 
Cappellari & Copin (120031) adaptive binning algorithm with the 
weighted Voronoi tessellation by Diehl & Statler C2006l see 
Rossetti & Molendi 120101 ). The temperature map is shown in 
Fig. |23] The formal errors given in output from the procedure 
are of the order of 15-20% but the measures at the map periph- 
ery, more affected by the background, should be considered with 
more caution. 
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Fig. 23. Temperature map (keV) with, superimposed, the contour 
levels of the XMM image. Coordinates are right ascension and 
dechnation (J2000). 



5. Discussion and conclusions 

5.1. Cluster structure 

The high values of the velocity dispersion cry = 1340^54^ 
km s"' and X-ray temperature kT = (6.38 + 0.25) keV are com- 
parable to the values found for hot, massive clusters (Mushotzky 
& Scharf [19971 Girardi & Mezzetti l200ll l. However, our esti- 
mates of (Tv and Tx are 3cr-inconsistent each other under the 
assumption of the equipartition of energy density between ICM 
and galaxies. In fact, we obtain ySspec - 1 ■71^0 2^ to be compared 
withySspec = 1, where ySspec = crl / (kT / fim^) v/ith fi - 0.58 the 
mean molecular weight and irip the proton mass (see also Fig.|6]l. 
Evidence that A2254 is well far from the dynamical equilibrium 
comes from both optical and X-ray analyses as discussed in the 
following. 

From the optical point of view, we find a strong evidence of 
bimodality from our 2D analysis of the photometric sample with 
the E and W-peaks far by -0.4-0.5 /i^,,' Mpc. The (modified) DS- 
test also shows evidence of substructure at a very high level (at 
the 99.99% c.L). Our 3D analysis of the region of the 2D E-peak, 
detects a high velocity subcluster which explains the high veloc- 
ity dispersion in that region and the strong increase of velocity 
dispersion profile at ~ 0.4-0.5 hj^ Mpc from the cluster center 
(see Fig.|6]l. 

However, a simple bimodal structure is not completely sat- 
isfactory when looking at the whole observational scenario. In 
fact, the W-peak is not well centered onto the BCG and both E- 
and W-peaks seem to have, at their turn, a bimodal or a N-S elon- 
gated structure (see Figs. [12] and [T3T l. The E-subcluster recov- 



ered through the 3D analysis is a minor structure and contains 
no dominant galaxy: it seems improbable that it causes the pres- 
ence of the 2D E-peak which is almost as rich as the W-peak (see 
Table[2]i. Indeed, more refined 3D analyses (e.g., the H-tree and 
the 3D-DEDICA methods) show evidence for a more complex 
structure. In particular, the H-tree method detects a low-velocity 
eastern subcluster which can explain even better the important 
2D E-peak and the strong increase of the velocity dispersion 
profile at ~ 0.4-0.5 /i^^q Mpc from the cluster center The low- 
velocity eastern subcluster might also explain the results of the 
weighted gap analysis (see Fig. [7]i. The 3D-DEDICA analysis 
reveals five subclusters, two of which in the western region: they 
might explain the presence of a bimodal 2D W-peak. 

We stress that our identification of subclusters are always 
confirmed by the 3D-KMM analysis. In principle, one might 
use the significance (as computed according to the likelihood 
ratio test statistics), to select the best cluster partition (e.g.. Bird 
1 19941 1. However, we note that in A2254, the high velocity E- 
group is the only subcluster well contained in the region sam- 
pled by our redshift data, while other subclusters have a large 
size and/or are at the boundary of this region. Thus, although the 
use of the 3D-KMM approach seems useful to perform a homo- 
geneous check of the subclusters detected using other disparate 
methods, we do not rely on the 3D-KMM results to fix the best 
partition of subclusters and the relevant kinematical properties. 
To obtain more stringent conclusions on the structure of A2254 
we stress the need of more galaxy redshifts and a more extended 
sampled region. 

As for the X-ray data, all the methods of the morphological 
analysis provide evidence of a dinamically disturbed cluster: i) 
the power ratio analysis, with its value of Pj'^/Pq'^ - 1.85 + 0.47 
together with the radio power at 1.4 GHz of Pi.4ghz = 2.9 x 
10^"^ h-^ W Hz ' (Govoni et al. 1200 Ibl l places A2554 in the re- 
gion of well known radio halo mergers, such as Abell 2256 and 
Abell 665, on the Pi/Po - Pi.4ghz relation of Buote (120011 see 
their Fig. 1-left, where Pi.4ghz = 2.1 x 10^"* W Hz"' in their cos- 
mology); ii) the combination of the three morphology indicators 
(as done by Cassano et al. I2010bl l in the (w, c), {P-ilPa,c), and 
{Pi/Po,c) planes designates A2254 as a morphologically dis- 
tubed cluster among the other clusters with radio halo in the sam- 
ple of that paper. In fact, A2254 shows well higher values of w 
and P3/P0 and a lower value of c with respect to the median dis- 
criminating values between disturbed and relaxed systems, i.e. 
w = 0.012 X 500 h-^ kpc, c = 0.2, and P3/P0 = 1.2 x 10"^ iii) 
The value of w = (0.037 + 0.001) x 500 /i^d kpc gives a sig- 
nificant indication of a non relaxed system. It is an intermediate 
value among those measured for 45 ROSAT-observed clusters 
by O'Hara et al. (2006). According to the simulations by Poole 
et al. (120061 their Fig. 11) this value of w is typical of ongoing 
cluster mergers. 

The active dynamical state of A2254 is also confirmed from 
the thermodinamic point of view, given the absence of a temper- 
ature drop in the core and the classification of A2254 as an high 
entropy core according to Leccardi et al. (120101 1 who fixed the 
separation value between medium and high entropy core clus- 
ter at £=0.64 (£= 0.833 + 0.097 in A2254). In fact, observa- 
tions seem to favor cool core destruction through cluster mergers 
(Allen et al .120011 Sanderson et al. 2006) and the vast majority of 
merging systems hosts high entropy cores according to Leccardi 
et al. ('2010), who analyzed 59 clusters spanning a large range of 
dynamical states. The temperature map, limited by the quality of 
the data to the inner regions, reveal a hot core with no significant 
temperature variations, consistent with an on-going merger 
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5.2. Optical - X-ray - radio data comparison 

At large scale, when comparing data from different wavelengths, 
the main feature is the correlation between the elongation of the 
optical structure, with the two galaxy clumps far by ~0.5 /i^q 
Mpc along the E-W direction, and the E-W elongation of the 
X-ray emission isophotes (see Fig.[T]l. The elongation of X-ray 
isophotes in the merging direction is expected from simulations 
(e.g., Roettiger et al. 1996 ). At small scale, we note that the X- 
ray emission is peaked at the BCG galaxy. Instead, the E-peak 
does not coincide with a secondary X-ray peak. This discrepancy 
between the distributions of the ICM (collisional) and galaxies 
(non collisional component) is a strong point in favor of an on- 
going cluster merger 

Our interpretation of the above correlations or no- 
correlations in the context of A2254 scenario is the following. 
The almost "triangular" shape of external X-ray isophotes in- 
dicates that the high velocity E-group has just passed through 
the main system going from west to east, as suggested by the 
comparison of Fig. |20]with, e.g., the second or third panels - or 
maybe the fourth panel considering a possible projection effect 
- in Figs. 4 and 5 of Poole et al. (2006 - top line, where the sub- 
cluster comes from the left). This scenario agrees with the fact 
that we observe the high-velocity E-group peak of galaxies being 
more advanced in the merging direction while the ICM compo- 
nent is slowed down and left behind (see e.g., the bullet cluster 
by Markevitch et al. 120021 1. At difference with the bullet cluster, 
A2254 does not reveals a secondary X-ray peak: this might be 
due to a more advanced phase of the merger or, in agreement 
with our mass estimates (see below), with the fact that the sub- 
cluster is much smaller and thus completely stripped by its gas 
content. 

However, a simple bimodal merger is not completely satis- 
factory. In fact, the X-ray emission in the internal cluster region, 
i.e. in the region of the W-peak, also reveals an elongation toward 
the North, perpendicular to the merging direction (see Figs. [T9l 
andl20t. This could be connected to a particular phase of the the 
E-W merger resembling the high emission structure perpendicu- 
lar to the merger at the time of the apocenter (but not in its sym- 
metric aspect, see Fig. 4 of Poole et al. 2006). Alternatively, it 
could be connected with another accretion phenomenon - along 
the N-S direction. In fact, optical 2D analyses of deep Subaru 
data (see Figs. [12] and Figs. [T3T l reveal a N-S double structure 
in the cluster center: one enclosing the BCG and one slightly at 
north. The 3D DEDICA analysis finds a northern and a southern 
galaxy clumps in the western part of the cluster, too. 

As for its radio morphology, A2254 is the most irregular in 
its radio morphology in the sample of Govoni et al. (20011?). 
The halo is not centered on the region of the maximum X-ray 
emission (see its upside-down "U" shape in Figs. [Tl and |20]|. re- 
sembling in some respects the behavior of radio relics. As for 
the comparison between radio and other data, the main feature is 
its bad correlation, already at large scales, with X-ray emission, 
and with the merging direction suggested by the two optical sub- 
clusters. In fact, here the halo does not permeate the cluster vol- 
ume in a similar way to the X-ray emitting gas of the intracluster 
medium (see Fig.[TJ as generally shown in other clusters not only 
by their large scale morphology but also by the point-to-point 
correlation between radio and X-ray brightness (see the case of 
the four clusters analyzed by Govoni et al. 12001 al l. Moreover, 
halos are round or elongated in the direction of the merger (the 
bullet cluster by Mai'kevitch et al. 20 02, Ab ell 520 by Girai'di et 
al.|2008i Abell 754 by Macario et al. 1201 11 1: this is not the case 
of the A2254 halo. 



Our interpretation is that the peculiar radio morphology of 
A2254 is likely due to the very complex substructure of A2254. 
In particular, the two N-S elongated branches of the radio emis- 
sion might be connected with the minor N-S mergers suggested 
by our optical analysis. However, note that, although the "U" 
shape morphology of the radio halo (see Fig. [T) recalls the "H' 
shape of the contours in the galaxy distribution in Fig. [12] the 
two features are shifted one respect to the other along the E-W 
direction. 

Finally, we stress the absence of a connection between bright 
galaxies and subclusters, while this is often found in complex 
systems in 2D (e.g ., Gell er & Beers 1982), in 3D (e.g., Abell 
520, Girardi et al. 120081 ). and also when a subcluster is par- 
tially destroyed (see the beautiful example of the tidal debris 
in A2744 by Owers et al. 120111 . In A2254, neither the high- 
velocity E-peak, the best detected through our optical analysis, 
hosts a bright, dominant galaxy. The feature which is worthy of 
mention, in spite of the poor statistics, is that BCG2, BCG3, and 
BCG4 are all located on the eastern side of the cluster, more 
external than respect to the E-subcluster However, their optical 
spectra do not show features - e.g. emission lines or a strong H^ 
absorption - typical of a possible starburst or poststarburst activ- 
ity triggered by the shocked gas during a cluster-cluster merger 
(Bekki et al. 120101 ). thus there is no clear connection with the 
ongoing merger. 

5.3. A simple model for mass and merging scenario 

Following the above discussion, A2254 can be described in first 
approximation by a main system and an eastern high velocity 
group. Looking at Table [3] we can estimate cry ~ 1000-1200 
km s"' and cry = 200-500 km s"' for the main and the E-group, 
respectively. 

Making the usual assumptions for the main system and the 
subcluster (cluster sphericity, dynamical equilibrium, coinci- 
dence of the galaxy and mass distributions), we can compute 
virial global quantities. We follow the prescriptions of Girardi & 
Mezzetti (2001 , see also Girardi et al. 1998 ) and compute /?vii - 
an estimate for R20Q - and the mass contained within this radius. 
In particular, we assume for the radius of the quasi-virialized 
region 7?vir = 0.17 x o-y/H{z) /i^,} Mpc(see Eq. 1 of Girardi & 
Mezz etti |200T] with the scaling of H(z) of Eq. 8 of Carlberg 
et al. [1997] for R20Q). For the mass we use M = My - S PT = 
in 12 ■ crlRpv/G - SPT (Eq. 3 of Girardi & Mezzetti |200TT ). 
with Rpv computed using the full procedure (see also Eq. 13 of 
Girardi et al. 1998 with A = /?vii for a close approximation) and 
the surface pressure term correction SPT = 0.2 x My. In prac- 
tice, both /?vir and M can be computed on the basis of the esti- 
mated velocity dispersion with the usual scaling-laws /?vii cry 
and M(< /?vii-) ^ o-y- 

We compute M(< /?vii- = 2.2-2.7 /i^^Mpc) = 1.5-2.5 
xlO'^ h-^ Mo and M(< R^i, = 0.4-1.2 /i^o'Mpc) = 0.1-1.8 
xlO'4 /i7,J Mq for the main and the E-group respectively, leading 
to a small mass ratio for the merger (< (1 : 10)). We compute a 
mass Msys(< R = 2.2-2.7 h-^Mpc) = 1.5-2.9 xlO'^ h'^ Mq for 
the whole system, where we assume that M cc Rto roughly ex- 
trapolate the mass of the E-group. According to the virial scal- 
ing laws, the above mass range leads to cry ~ (1000-1250) 
km s ' for the virial velocity dispersion, thus making acceptable 
the ySspec = 1 model. 

For comparison, we would obtain Msys,i(< /?vii,i = 
3.0 /i^^Mpc) = (3.5 + 1.0) xlO'^ h-^ Mq when assum- 
ing A2254 as a dynamically relaxed cluster and Msys,5 ~ 4 
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xlO'^ h.-^ Mq when assuming A2254 is the very complex clus- 
ter formed by the five DEDICA subclusters. 

We use X-ray data to compute an alternative mass estimate. 
Given the no significant departures from the isothermality of 
the temperature profile out to ~ 0.9 /i^q Mpc (see Fig. l22l l. we 
derive a mass estimate using the isothermal beta model where 
the total mass profile can be expressed by a simple analyti- 
cal formula (e.g., eq. 6 of Henry et al. 1 19931 with y = 1). 
Using the value of ^rouT and the fitted parameters of the /3 
model to the surface brightness (Sect. 14. Il l, we obtain a mass 
Mx(< R = 0.9 /i^fjMpc) = 0.48 xlO'^ /i^,,' with nominal er- 
rors less than the 5%. We rescale the optical mass Msy, toR = 0.9 
h^jQ Mpc assuming that the system is described by a King-like 
mass distribution with a very small core radius (Girardi et al. 
11998 ) or, alternatively, a NFW profile where the mass-dependent 
concentration parameter c is taken from Navarro et al. (1997]l 
and rescaled by the factor I + z (Bullock et al. 1200 II Dolag et al. 
I2004D . here c ~ 4. We obtain M,y,{< R = 0.9 /z^oMpc) = (0.6- 
1.0) xlO'^ /z^o Mq to be compared with the above Mx(< R - 
0.9 /i^fjMpc) ~ 0.5 xlO'^ Mq. Thus, when opportunely cor- 
recting the optical mass for the presence of substructure, we do 
not find a large discrepancy between optical and X-ray mass esti- 
mates in agreement with thatySspec = 1 is acceptable (see above). 

Taking the above face on values, one might estimate that the 
non-thermal pressure in A2254 ranges from a value of 10% to 
100% of the gas thermal pressure. However, note that the above 
comparison of the two mass estimates should be considered at 
best a tentative since the multiclump structure of the cluster 
should be better known. Indeed, previous efforts on making a 
comparison to put a limit on non-thermal pressure have targeted 
systems which can be assumed in equilibrium at least in the non- 
collisional component (e.g., a sample of elliptical galaxies by 
Churazov et al. l20T0] l. 

We can use our results to examine A2254 with respect to the 
observed scaling relations among Piaghz, the halo radio power, 
and Mh, the total cluster mass contained within the size of the 
radio halo, Rh (Cassano et al. ,2007) . We rescale the mass Msy, 
to Rn = 0.4 /i7(5 Mpc (Cassano et al. I2007I I and obtain Msys(< 
Rh = 0.4 /!7,!Mpc) = (0.2-0.5) xlO^^ /i^d Mq. Thus A2254 is 
well in agreement with the relation fitted for other clusters with 
radio halos (see Fig. 9 of Cassano et al. l2007l l. 

We also attempt to analyze the cluster merger between the 
main system and the high velocity E-group. Even if the sub- 
structure of A2254 is more complex, the elongation of the X- 
ray surface brightness suggests that this merger is the most re- 
cent one. We apply the two-body model (Beers et al. 119821 
Thompson 1 19821 1 following the methodology outlined for other 
DARC clusters, e.g. Abell 520 (Girardi et al. 2008) and Abell 
2345 (Boschin et al. 120101 ). The values of the relevant parame- 
ters for the two-subclusters system are: the above Mjys - 1.5- 
2.9 xlO'^ /j^g Mq, the relative LOS velocity in the rest-frame, 
AVrf.LOS ~ 3000 km s"\ and the projected linear distance be- 
tween the two clumps, D ~ 0.5 /j^g Mpc, taken as the distance 
between the BCG and the 2D E-peak. The comparison between 
the observed shape of the X-ray surface brightness isophotes and 
the simulated ones (Poole et al. 2006) suggests we can assume 
to see the merger after the core crossing and not too much af- 
ter the apocenter For comparison with other clusters hosting 
radio halos, the time t elapsed from the core crossing is gen- 
erally few fractions of Gyr (e.g., Markevitch et al. 2002; Girardi 
et al. 120081) . Figure l24l compares the bimodal-model solutions 
as a function of a, where a is the projection angle between the 
plane of the sky and the line connecting the centers of the two 
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Fig. 24. System mass vs. projection angle for bound and un- 
bound solutions (thick solid and thick dashed curves, respec- 
tively) of the two-body model applied to the high velocity E- 
group and the main cluster Thin/black and thick/red lines refer 
to f = 0.1 Gyr and f = 0.5 Gyr, respectively. Labels Bla and Bib 
indicate the bound and incoming, i.e., collapsing solutions (solid 
curve). Label BO and UO indicate the bound outgoing, i.e., ex- 
panding solutions and unbound outgoing solutions (solid curve 
going on in the dotted curve, respectively). The horizontal lines 
give the range of observational values of the mass system. The 
thin dashed curve separates bound and unbound regions accord- 
ing to the Newtonian criterion (above and below the thin dashed 
curve, respectively). 



clumps, with the observed mass of the system Msys. We con- 
sider two different cases (f = 0.1 Gyr and t - 0.5 Gyr). In both 
cases we find a bound outgoing solution (BO) with a ~ 45°, and 
a ~ 75°, respectively. These intermediate angles are someway 
in agreement with the fact that we can detect the E-peak in the 
3D analysis, which is the most efficient for intermediate angles 
(Pinkney et al. 1996). These angles lead to a deprojected relative 
velocity AVrf ~ 3100^200 km s"' of the subcluster with re- 
spect to the main system, comparable to that of the bullet cluster 
(Markevitch et al. 2002 ) or of Abell 2744 (Boschin et al. 2006)). 

In conclusion, A2254, for its large mass and unrelaxed dy- 
namical status, fits well among typical clusters with radio halos 
described in the literature. Present data agree with the scenario of 
a recent merger with a minor subcluster, but we have evidence 
for a more complex dynamical history which should be better 
investigated to explain the peculiar radio morphology. 
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Fig. 1. Subaru V-band image of the cluster A2254 (North at the top and East to the left) with, superimposed, the contour levels of 
the XMM archival image (thick/blue contours; photons in the energy range 0.5-2 keV) and the contour levels of a VLA radio image 
at 1.4 GHz (thin/red contours, see Govoni et al. 2001b). Labels indicate galaxies cited in the text. The two (green) circles indicate 
the positions of the E- and W-peaks in the galaxy distribution from the 2D-DEDICA analysis (Subaru sample, with /' < 20.5 mag, 
see Tab. IS . 
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Fig. 2. INT r'-band image of the cluster A2254 (North at the top and East to the left). Circles and squares indicate cluster members 
and non-members, respectively (see Table[TJ. Labels indicate cluster galaxies cited in the text. 
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Fig. 3. Subaru /'-band image of the very central cluster region. 
It highlights the BCG and its bright spiral companion BCGc 
(North at the top and East to the left). 




Fig. 4. Multiwavelength picture of the region around RG galaxy. 
The MlPS-Spitzer image (green color) is superimposed to the 
Subaru /'-band image. The arrow indicates the faint extended 
object close to RG. The black circle indicate the position of the 
XMM point source. Contour levels of the radio image are shown 
in blue (Govoni et al. 1200 Ibl l. North at the top and East to the 
left. 



